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1 Introduction 

White dwarf binaries are thought to be the most common binaries in the Uni- 
verse, and in our Galaxy their number is estimated to be as high as 10 8 . In 
addition most stars are known to be part of binary systems, roughly half of 
which have orbital periods short enough that the evolution of the two stars 
is strongly influenced by the presence of a companion. Furthermore, it has 
become clear from observed close binaries, that a large fraction of binaries 
that interacted in the past must have lost considerable amounts of angular 
momentum, thus forming compact binaries, with compact stellar components. 
The details of the evolution leading to the loss of angular momentum are un- 
certain, but generally this is interpreted in the framework of the so called 
"common-envelope evolution" : the picture that in a mass-transfer phase be- 
tween a giant and a more compact companion the companion quickly ends 
up inside the giant's envelope, after which frictional processes slow down the 
companion and the core of the giant, causing the "common envelope" to be 
expelled, as well as the orbital separation to shrink dramatically [55] . 

Among the most compact binaries know, often called ultra-compact or 
ultra-short binaries, are those hosting two white dwarfs and classified into two 
types: detached binaries, in which the two components are relatively widely 
separated and interacting binaries, in which mass is transferred from one 
component to the other. In the latter class a white dwarf is accreting from a 
white dwarf like object (we often refer to them as AM CVn systems, after the 
prototype of the class, the variable star AM CVn; [56, 28]). 

In the past many authors have emphasised the importance of studying 
white dwarfs in DDBs. In fact, the study of ultra-short white dwarf binaries is 
relevant to some important astrophysical questions which have been outlined 
by several author. Recently, [32] listed the following ones: 

• Binary evolution Double white dwarfs are excellent tests of binary evo- 
lution. In particular the orbital shrinkage during the common-envelope 
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Fig. 1. Period versus total mass of double white dwarfs. The points and arrows are 
observed systems [29], the grey shade a model for the Galactic population. Systems 
to the left of the dashed line will merge within a Hubble time, systems above the 
dotted line have a combined mass above the Chandrasekhar mass. The top left 
corner shows the region of possible type la supernova progenitors, where the grey 
shade has been darkened for better visibility (adapted from [32]). 



phase can be tested using double white dwarfs. The reason is that for gi- 
ants there is a direct relation between the mass of the core (which becomes 
a white dwarf and so its mass is still measurable today) and the radius of 
the giant. The latter carries information about the (minimal) separation 
between the two components in the binary before the common envelope, 
while the separation after the common envelope can be estimated from 
the current orbital period. This enables a detailed reconstruction of the 
evolution leading from a binary consisting of two main sequence stars to 
a close double white dwarf [26] . The interesting conclusion of this exercise 
is that the standard schematic description of the common envelope - in 
which the envelope is expelled at the expense of the orbital energy - can- 
not be correct. An alternative scheme, based on the angular momentum, 
for the moment seems to be able to explain all the observations [30] . 
• Type la supernovae Type la supernovae have peak brightnesses that are 
well correlated with the shape of their light curve [35] , making them ideal 
standard candles to determine distances. The measurement of the appar- 
ent brightness of far away supernovae as a function of redshift has led to 
the conclusion that the expansion of the universe is accelerating [34, 45]. 
This depends on the assumption that these far-away (and thus old) super- 
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novae behave the same as their local cousins, which is a quite reasonable 
assumption. However, one of the problems is that we do not know what 
exactly explodes and why, so the likelihood of this assumption is difficult 
to assess [37]. One of the proposed models for the progenitors of type la 
supernovae are massive close double white dwarfs that will explode when 
the two stars merge [16]. In Fig. 1 the observed double white dwarfs are 
compared to a model for the Galactic population of double white dwarfs 
[27], in which the merger rate of massive double white dwarfs is similar to 
the type la supernova rate. The grey shade in the relevant corner of the 
diagram is enhanced for visibility. The discovery of at least one system in 
this box confirms the viability of this model (in terms of event rates) . 
Accretion physics The fact that in AM CVn systems the mass losing star 
is an evolved, hydrogen deficient star, gives rise to a unique astrophysical 
laboratory, in which accretion discs made of almost pure helium [22, 49, 
13, 46, 57]. This opens the possibility to test the behaviour of accretion 
discs of different chemical composition. 

Gravitational wave emission Untill recently the DDBs with two NSs were 
considered among the best sources to look for gravitational wave emission, 
mainly due to the relatively high chirp mass expected for these sources, 
In fact, simply inferring the strength of the gravitational wave amplitude 
expected for from [11] 
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where the frequency of the wave is given by / = 2/P or i,. It is evident that 
the strain signal h from DDBs hosting neutron stars is a factor 5-20 higher 
than in the case of DDBs with white dwarfs as far as the orbital period is 
larger than approximativcly 10-20 minutes. In recent years, AM CVns have 
received great attention as they represent a large population of guaranteed 
sources for the forthcoming Laser Interferometer Space Antenna [31, 51]. 
Double WD binaries enter the LISA observational window (0.1 100 
mHz) at an orbital period ~ 5 hrs and, as they evolve secularly through 
GW emission, they cross the whole LISA band. They are expected to be so 
numerous (~ 10 3 4- 10 4 expected), close on average, and luminous in GWs 
as to create a stochastic foreground that dominates the LISA observational 
window up to « 3 mHz [51]. Detailed knowledge of the characteristics of 
their background signal would thus be needed to model it and study weaker 
background GW signals of cosmological origin. 
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Fig. 2. Expected signals of ultra-compact binaries, the ones with error bars from 
(adapted from [46, 32]. 

A relatively small number of ultracompact DDBs systems is presently 
known. According to [42] there exist 17 confirmed objects with orbital pe- 
riods in the 10 -j- 70 min in which a hydrogen-deficient mass donor, either a 
semi-degenerate star or a WD itself, is present. These are called AM CVn 
systems and are roughly characterized by optical emission modulated at the 
orbital period, X-ray emission showing no evidence for a significant modula- 
tion (from which a moderately magnetic primary is suggested, [43]) and, in 
the few cases where timing analyses could be carried out, orbital spin-down 
consistent with GW emission-driven mass transfer. 

In addition there exist two peculiar objects, sharing a number of obser- 
vational properties that partially match those of the "standard" AM CVn's. 
They are RX J0806.3+1527 and RX J1914.4+2456, whose X-ray emission is 
~ 100% pulsed, with on-phase and off-phase of approximately equal duration. 
The single modulations found in their lightcurves, both in the optical and in 
X-rays, correspond to periods of, respectively, 321.5 and 569 s ([21, 52]) and 
were first interpreted as orbital periods. If so, these two objects are the binary 
systems with the shortest orbital period known and could belong to the AM 
CVn class. However, in addition to peculiar emission properties with respect 
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to other AM CVn's, timing analyses carried out by the above cited authors 
demonstrate that, in this interpretation, these two objects have shrinking or- 
bits. This is contrary to what expected in mass transferring double white 
dwarf systems (including AM CVn's systems) and suggests the possibility 
that the binary is detached, with the orbit shrinking because of GW emission. 
The electromagnetic emission would have in turn to be caused by some other 
kind of interaction. 

Nonetheless, there are a number of alternative models to account for the 
observed properties, all of them based upon binary systems. The intermediate 
polar (IP) model ([25, 17, 33]) is the only one in which the pulsation periods 
are not assumed to be orbital. In this model, the pulsations are likely due 
to the spin of a white dwarf accreting from non-degenerate secondary star. 
Moreover, due to geometrical constraints the orbital period is not expected 
to be detectable. The other two models assume a double white dwarf binaries 
in which the pulsation period is the orbital period. Each of them invoke a 
semi-detached, accreting double white dwarfs: one is magnetic, the double 
degenerate polar model ([6, 39, 40, 18]), while the other is non-magnetic, the 
direct impact model ([27, 24, 39]), in which, due to the compact dimensions 
of these systems, the mass transfer streams is forced to hit directly onto the 
accreting white dwarfs rather than to form an accretion disk . 

After a brief presentation of the two X-ray selected double degenerate bi- 
nary systems, we discuss the main scenario of this type, the Unipolar Inductor 
Model (UIM) introduced by [59] and further developed by [7, 8], and compare 
its predictions with the salient observed properties of these two sources. 

1.1 RXJ0806. 3+1527 

RXJ0806.3+1527 was discovered in 1990 with the ROSAT satellite during 
the All-Sky Survey (RASS; [3]). However, it was only in 1999 that a periodic 
signal at 321s was detected in its soft X-ray flux with the ROSAT HRI ([17, 
4]). Subsequent deeper optical studies allowed to unambiguously identify the 
optical counterpart of RXJ0806. 3+1527, a blue V = 21.1 (B = 20.7) star 
([18, 19]). B, V and R time-resolved photometry revealed the presence of a 
~ 15% modulation at the ~ 321s X-ray period ([19, 39]. 

The VLT spectral study revealed a blue continuum with no intrinsic ab- 
sorption lines [19] . Broad (FWHM ~ 1500 km s _1 ), low equivalent width 
(EW ~ —2 -. — 6 A) emission lines from the He II Pickering series (plus ad- 
ditional emission lines likely associated with He I, C III, N III, etc.; for a 
different interpretation see [44]) were instead detected [19]. These findings, 
together with the period stability and absence of any additional modulation 
in the 1 min-5 hr period range, were interpreted in terms of a double degen- 
erate He- rich binary (a subset of the AM CVn class; see [56]) with an orbital 
period of 321s, the shortest ever recorded. Moreover, RX J0806. 3+1527 was 
noticed to have optical/X-ray properties similar to those of RX J1914. 4+2456, 
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Table 1. Overview of observational properties of AM CVn stars (adapted from [28]) 
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a orb = orbital, sh = superhump, periods from ww03, see references therein, 
(p)/(s)/(X) for photometric, spectroscopic, X-ray period. 
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c A = ASCA, C = Chandra, R = ROSAT, Rx = RXTE, X = XMM-Newton kns+04 
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a 569 s modulated soft X-ray source proposed as a double degenerate system 
([6, 38, 40]). 

In the past years the detection of spin-up was reported, at a rate of 
~6.2xl0~ n s s _1 , for the 321s orbital modulation, based on optical data 
taken from the Nordic Optical Telescope (NOT) and the VLT archive, and by 
using incoherent timing techniques [14, 15]. Similar results were reported also 
for the X-ray data (ROSAT and Chandra; [53]) of RX J0806.3+1527 spanning 
over 10 years of uncoherent observations and based on the NOT results [14]. 

A Tclescopio Nazionale Galileo (TNG) long-term project (started on 2000) 
devoted to the study of the long-term timing properties of RX J0806. 3+1527 
found a slightly energy-dependent pulse shape with the pulsed fraction in- 
creasing toward longer wavelengths, from ~12% in the B-band to nearly 14% 
in the I- band (see lower right panel of Figure 5; [21]). An additional variabil- 
ity, at a level of 4% of the optical pulse shape as a function of time (see upper 
right panel of Figure 5 right) was detected. The first coherent timing solution 
was also inferred for this source, firmly assessing that the source was spinning- 
up: P=321.53033(2)s, and P=-3.67(l) x lO" 11 s s" 1 (90% uncertainties are re- 
ported; [21]). Reference [54] obtained independently a phase-coherent timing 
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Fig. 3. VLT FORSl medium (6A; 3900-6000A) and low (30A; above 6000A) reso- 
lution spectra obtained for the optical counterpart of RX J0806. 3+1527. Numerous 
faint emission lines of Hel and Hell (blended with H) are labeled (adapted form 
[19]). 



solutions for the orbital period of this source over a similar baseline, that is 
fully consistent with that of [21]. See [2] for a similar coherent timing solution 
also including the covariance terms of the fitted parameters. 

The relatively high accuracy obtained for the optical phase coherent P-P 
solution (in the January 2001 - May 2004 interval) was used to extend its valid- 
ity backward to the ROSAT observations without loosing the phase coherency, 
i.e. only one possible period cycle consistent with our P-P solution. The best 
X-ray phase coherent solution is P=321.53038(2)s, P=-3.661(5) x lO" 11 s s" 1 
(for more details see [21]). Figure 5 (left panel) shows the optical (2001-2004) 
and X-ray (1994-2002) light curves folded by using the above reported P- 
P coherent solution, confirming the amazing stability of the X-ray/optical 
anti-correlation first noted by ([20]; see inset of left panel of Figure 5). 

On 2001, a Chandra observation of RXJ0806. 3+1527 carried out in simul- 
taneity with time resolved optical observation at the VLT, allowed for the first 
time to study the details of the X-ray emission and the phase-shift between 
X-rays and optical band. The X-ray spectrum is consistent with a occulting, 
as a function of modulation phase, black body with a temperature of ^60 eV 
[20]. A 0.5 phase-shift was reported for the X-rays and the optical band [20]. 
More recently, a 0.2 phase-shift was reported by analysing the whole historical 
X-ray and optical dataset: this latter result is considered the correct one [2]. 

On 2002 November 1 st a second deep X-ray observation was obtained with 
the XMM-Newton instrumentations for about 26000 s, providing an increased 
spectral accuracy (see eft panel of Figure 7). The XMM-Newton data show a 
lower value of the absorption column, a relatively constant black body temper- 
ature, a smaller black body size, and, correspondingly, a slightly lower flux. 
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Fig. 4. Left panel: Results of the phase fitting technique used to infer the P-P 
coherent solution for RXJ0806. 3+1527: the linear term (P component) has been 
corrected, while the quadratic term (the P component) has been kept for clarity. 
The best P solution inferred for the optical band is marked by the solid fit line. 
Right panel: 2001-2004 optical flux measurements at fdifferent wavelengths. 



All these differences may be ascribed to the pile-up effect in the Chandra 
data, even though we can not completely rule out the presence of real spec- 
tral variations as a function of time. In any case we note that this result is in 
agreement with the idea of a self-eclipsing (due only to a geometrical effect) 
small, hot and X-ray emitting region on the primary star. Timing analysis did 
not show any additional significant signal at periods longer or shorter than 
321.5s, (in the 5hr-200ms interval). By using the XMM-Newton OM a first 
look at the source in the UV band (see right panel of Figure 7) was obtained 
confirming the presence of the blackbody component inferred from IR/optical 
bands. 

Reference [20] measured an on-phase X-ray luminosity (in the range 0.1- 
2.5 keV) L x = 8 x 10 31 (d/200 pc) 2 erg s" 1 for this source. These authors 
suggested that the bolometric luminosity might even be dominated by the 
(unseen) value of the UV flux, and reach values up to 5-6 times higher. The 
optical flux is only ~ 15% pulsed, indicating that most of it might not be 
associated to the same mechanism producing the pulsed X-ray emission (pos- 
sibly the cooling luminosity of the WD plays a role). Given these uncertain- 
ties and, mainly, the uncertainty in the distance to the source, a luminosity 
W ~ 10 32 (d/200 pc) 2 erg s _1 will be assumed as a reference value. 
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Fig. 5. Left Panel: The 1994-2002 phase coherently connected X-ray folded light 
curves (filled squares; 100% pulsed fraction) of RX J0806. 3+1527, together with 
the VLT-TNG 2001-2004 phase connected folded optical light curves (filled circles). 
Two orbital cycles are reported for clarity. A nearly anti-correlation was found. 
Right panels: Analysis of the phase variations induced by pulse shape changes in 
the optical band (upper panel), and the pulsed fraction as a function of optical 
wavelengths (lower panel). 



1.2 RX J1914. 4+2456 

The luminosity and distance of this source have been subject to much de- 
bate over the last years. Reference [59] refer to earlier ASCA measurements 
that, for a distance of 200-500 pc, corresponded to a luminosity in the range 
(4 x 10 33 + 2.5 x 10 34 ) erg s _1 . Reference [41], based on more recent XMM- 
Newton observations and a standard blackbody fit to the X-ray spectrum, 
derived an X-ray luminosity of ~ 10 35 c^ pc erg s _1 , where d kpc is the distance 
in kpc. The larger distance of ~ 1 kpc was based on a work by [50]. Still more 
recently, [42] find that an optically thin thermal emission spectrum, with an 
edge at 0.83 keV attributed to O VIII, gives a significantly better fit to the 
data than a blackbody model. The optically thin thermal plasma model im- 
plies a much lower bolometric luminosity of L bol ~ 10 33 d^ pc erg s _1 . 
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Phase 

Fig. 6. The results of the XMM-Newton phase-resolved spectroscopy (PRS) analysis 
for the absorbed blackbody spectral parameters: absorption, blackbody temperature, 
blackbody radius (assuming a distance of 500 pc), and absorbed (triangles) and 
unabsorbed (asterisks) flux. Superposed is the folded X-ray light curve. 

Reference [42] also note that the determination of a 1 kpc distance is not 
free of uncertainties and that a minimum distance of ~ 200 pc might still be 
possible: the latter leads to a minimum luminosity of ~ 3 x 10 31 erg s _1 . 
Given these large discrepancies, interpretation of this source's properties re- 
mains ambiguous and dependent on assumptions. In the following, we refer 
to the more recent assessment by [42] of a luminosity L = 10 33 erg s _1 for a 
1 kpc distance. 

Reference [42] also find possible evidence, at least in a few observations, of 
two secondary peaks in power spectra. These are very close to (Av ~ 5 x 10~ 5 
Hz) and symmetrically distributed around the strongest peak at ~ 1.76 x 10~ 3 
Hz. References [42] and [10] discuss the implications of this possible finding. 
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Fig. 7. Broad-band energy spectrum of RX J0806. 3+1527 as inferred from the 
Chandra, XMM-Newton, VLT and TNG measurements and EUVE upper limits. 
The dotted line represents one of the possible fitting blackbody models for the 
IR/optical/UV bands. 



2 The Unipolar Inductor Model 

The Unipolar Inductor Model (UIM) was originally proposed to explain the 
origin of bursts of decametric radiation received from Jupiter, whose prop- 
erties appear to be strongly influenced by the orbital location of Jupiter's 
satellite Io [12, 36] . 

The model relies on Jupiter's spin being different from the system orbital pe- 
riod (Io spin is tidally locked to the orbit). Jupiter has a surface magnetic field 
~ 10 G so that, given Io's good electrical conductivity (a), the satellite expe- 
riences an e.m.f. as it moves across the planet's field lines along the orbit. The 
e.m.f. accelerates free charges in the ambient medium, giving rise to a flow of 
current along the sides of the flux tube connecting the bodies. Flowing charges 
are accelerated to mildly relativistic energies and emit coherent cyclotron ra- 
diation through a loss cone instability (cfr. [58] and references therein): this 
is the basic framework in which Jupiter decametric radiation and its modu- 
lation by Io's position are explained. Among the several confirmations of the 
UIM in this system, HST UV observations revealed the localized emission 
on Jupiter's surface due to flowing particles hitting the planet's surface - the 
so-called Io's footprint [5]. In recent years, the complex interaction between 
Io-related free charges (forming the Io torus) and Jupiter's magnetosphere 
has been understood in much greater detail [47, 48]. Despite these significant 
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complications, the above scenario maintains its general validity, particularly 
in view of astrophysical applications. 

Reference [59] considered the UIM in the case of close white dwarf bina- 
ries. They assumed a moderately magnetized primary, whose spin is not syn- 
chronous with the orbit and a non-magnetic companion, whose spin is tidally 
locked. They particularly highlight the role of ohmic dissipation of currents 
flowing through the two WDs and show that this occurs essentially in the 
primary atmosphere. A small bundle of field lines leaving the primary sur- 
face thread the whole secondary. The orbital position of the latter is thus 
"mapped" to a small region onto the primary's surface; it is in this small re- 
gion that ohmic dissipation - and the associated heating - mainly takes place. 
The resulting geometry, illustrated in Fig. 8, naturally leads to mainly ther- 
mal, strongly pulsed X-ray emission, as the secondary moves along the orbit. 
The source of the X-ray emission is ultimately represented by the relative mo- 



Fig. 8. Electric coupling between the asynchronous, magnetic primary star and the 
non-magnetic secondary,, in the UIM (adapted from [59]). 

tion between primary spin and orbit, that powers the electric circuit. Because 
of resistive dissipation of currents, the relative motion is eventually expected 
to be cancelled. This in turn requires angular momentum to be redistributed 
between spin and orbit in order to synchronize them. The necessary torque is 
provided by the Lorentz force on cross-field currents within the two stars. 
Reference [59] derived synchronization timescales (t q ) ~ few 10 3 yrs for both 
RXJ1914.4+2456 and RX J0806.3+1527 , less than 1% of their orbital evo- 
lutionary timescales. This would imply a much larger Galactic population of 
such systems than predicted by population-synthesis models, a major difficulty 
of this version of the UIM. However, [7, 8] have shown that the electrically ac- 
tive phase is actually long-lived because perfect synchronism is never reached. 
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In a perfectly synchronous system the electric circuit would be turned off, 
while GWs would still cause orbital spin-up. Orbital motion and primary spin 
would thus go out of synchronism, which in turn would switch the circuit on. 
The synchronizing (magnetic coupling) and de-synchronizing (GWs) torques 
are thus expected to reach an equilibrium state at a sufficiently small degree 
of asynchronism. 

We discuss in detail how the model works and how the major observed prop- 
erties of RXJ0806.3+1527 and RX J1914.4+2456 can be interpreted in the 
UIM framework. We refer to [1] for a possible criticism of the model based on 
the shape of the pulsed profiles of the two sources. Finally, we refer to [9, 10], 
who have recently proposed alternative mass transfer models that can also 
account for long-lasting episodes of spin-up in Double White Dwarf systems. 



3 UIM in Double Degenerate Binaries 

According to [59], define the primary's asynchronism parameter a = oj\/ui , 
where uj\ and uj are the primary's spin and orbital frequencies. In a system 
with orbital separation a, the secondary star will move with the velocity v = 
a{bj — bJi) = [GM\(1 + q)} 1 ^ 3 w<!/ 3 (l — a) relative to field lines, where G is the 
gravitational constant, Mi the primary mass, q — M^/Mi the system mass- 
ratio. The electric field induced through the secondary is thus E = — -, 
with an associated e.m.f. <P — 2i? 2 -E, Ri being the secondary's radius and B 2 
the primary magnetic field at the secondary's location. The internal (Lorcntz) 
torque redistributes angular momentum between spin and orbit conserving 
their sum (see below) , while GW-emission causes a net loss of orbital angular 
momentum. Therefore, as long as the primary spin is not efficiently affected by 
other forces, i.e. tidal forces (cfr. App.A in [7]), it will lag behind the evolving 
orbital frequency, thus keeping electric coupling continuously active. 
Since most of the power dissipation occurs at the primary atmosphere (cfr. 
[59]), we slightly simplify our treatment assuming no dissipation at all at 
the secondary. In this case, the binary system is wholly analogous to the 
elementary circuit of Fig. 9. Given the e.m.f. (<£) across the secondary star 
and the system's effective resistivity TZ « (2cri? 2 ) _1 (a/i?i) 3 / 2 , the dissipation 
rate of electric current (W) in the primary atmosphere is: 

W = I 2 1l = I$ = kLjl 7 / 3 (l - a) 2 (4) 

where k = 2(fi 1 / l c) 2 aR\ /2 : ^/[GMi(l + q)} 11/6 is a constant of the system. 
The Lorentz torque (Nl) has the following properties: i) it acts with the 
same magnitude and opposite signs on the primary star and the orbit, Nl — 
N { L 1] = -N L orh) . Therefore; ii) N L conserves the total angular momentum in 
the system, transferring all that is extracted from one component to the other 
one; Hi) Nl is simply related to the energy dissipation rate: W — NLL0 o (l — a). 
From the above, the evolution equation for lo\ is: 
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Fig. 9. Sketch of the elementary circuit envisaged in the UIM. The secondary star 
acts as the battery, the primary star represents a resistance connected to the battery 
by conducting "wires" (field lines.) Inclusion of the effect of GWs corresponds to 
connecting the battery to a plug, so that it is recharged at some given rate. Once 
the battery initial energy reservoir is consumed, the bulb will be powered just by 
the energy fed through the plug. This corresponds to the "steady-state" solution. 
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The orbital angular momentum is L a = I u , so that the orbital evolution 
equation is: 



3{N e „ + N { L ° rb) )/I = -3(N^-N L )/I = 



E„- 
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where I = q(l + q) 1 / 3 G 2 / 3 Af^ 3 tj 4 ^ 3 is the orbital moment of inertia and 
N s „ = E s „/uj is the GW torque. 



3.1 Energetics of the electric circuit 

Let us focus on how energy is transferred and consumed by the electric circuit. 
We begin considering the rate of work done by Nl on the orbit 



E [ ° rh) = N^ujo 



and that done on the primary: 

E S pin = N^UJi = 



-NtM 



w 

~ a' 



1 — a 



-W 



-aE 



(orb) 



(7) 



(8) 



The sum E. 



E ( ° rb) — -W. Clearly, not all of the energy extracted from 
one component is transferred to the other one. The energy lost to ohmic 



spin 
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dissipation represents the energetic cost of spin-orbit coupling. 
The above formulae allow to draw some further conclusions concerning the 
relation between a and the energetics of the electrical circuit. When a > 
1, the circuit is powered at the expenses of the primary's spin energy. A 
fraction a -1 of this energy is transferred to the orbit, the rest being lost to 
ohmic dissipation. When a < 1, the circuit is powered at the expenses of the 
orbital energy and a fraction a of this energy is transferred to the primary 
spin. Therefore, the parameter a represents a measure of the energy transfer 
efficiency of spin-orbit coupling: the more asynchronous a system is, the less 
efficiently energy is transferred, most of it being dissipated as heat. 



3.2 Stationary state: General solution 

As long as the asynchronism parameter is sufficiently far from unity, its evo- 
lution will be essentially determined by the strength of the synchronizing 
(Lorentz) torque, the GW torque being of minor relevance. The evolution in 
this case depends on the initial values of a and u> , and on stellar parameters. 
This evolutionary phase drives a towards unity, i.e. spin and orbit are driven 
towards synchronism. It is in this regime that the GW torque becomes im- 
portant in determining the subsequent evolution of the system. 
Once the condition a = 1 is reached, indeed, GW emission drives a small 
angular momentum disequilibrium. The Lorentz torque is in turn switched on 
to transfer to the primary spin the amount of angular momentum required 
for it to keep up with the evolving orbital frequency. This translates to the 
requirement that Cj\ = uj . By use of expressions (5) and (6), it is found that 
this condition implies the following equilibrium value for a (we call it ctoo): 

°° ~ k W "/3 W 

This is greater than zero if the orbit is shrinking (6j > 0), which implies that 
ctoo < 1- For a widening orbit, on the other hand, > 1. However, this 
latter case does not correspond to a long-lived configuration. Indeed, define 
the electric energy reservoir as Ejjim = {\/2)I\(u)\ — u>%), which is negative 
when a < 1 and positive when a > 1. Substituting eq. (9) into this definition: 

Euim = -W, (10) 

If 

a — a oo > 1; energy is consumed at the rate W: the circuit will eventually 
switch off (aoo = !)• At later times, the case < 1 applies. 
If a — aoo < 1, condition (10) means that the battery recharges at the rate 
W at which electric currents dissipate energy: the electric energy reservoir is 
conserved as the binary evolves. 

The latter conclusion can be reversed (cfr. Fig. 9): in steady-state, the rate 
of energy dissipation (W) is fixed by the rate at which power is fed to the 
circuit by the plug (E UIM ). The latter is determined by GW emission and 
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the Lorentz torque and, therefore, by component masses, u and fi]_. 
Therefore the steady-state degree of asynchronism of a given binary system 
is uniquely determined, given oj . Since both oj and 6j evolve secularly, the 
equilibrium state will be "quasi-steady" , evolving secularly as well. 

3.3 Model application: equations of practical use 

We have discussed in previous sections the existence of an asymptotic regime 
in the evolution of binaries in the UIM framework. Given the definition of a x 
and W (eq. 9 and 4, respectively), we have: 

W = hUoLJo Q - ^ ■ (11) 

The quantity (1 — aoo) represents the actual degree of asynchronism only 
for those systems that had enough time to evolve towards steady-state, i.e 
with sufficiently short orbital period. In this case, the steady-state source 
luminosity can thus be written as: 

Woo = /iw w (l - aoo) (12) 

Therefore - under the assumption that a source is in steady-state - the quantity 
aoo can be determined from the measured values of W, uj , lu . Given its 
definition (eq. 9), this gives an estimate of k and, thus, fi\. 
The equation for the orbital evolution (6) provides a further relation between 
the three measured quantities, component masses and degree of asynchronism. 
This can be written as: 

£ gr + \i o u 2 {lu /u ) = — (13) 
A (1 — a) 

that becomes, inserting the appropriate expressions for E ST and I a : 

"^^-i^X + ^-O. (14, 
where X = M^' 3 q/(1 + q) 1 ^ 3 = M. 5 / 3 , M. being the system's chirp mass. 



4 RXJ0806.3+1527 

We assume here the values of and of the bolomctric luminosity reported 

in 1 and refer to See [8] for a complete discussion on how our conclusions 
depend on these assumptions. 

In Fig. 10 (see caption for further details), the dashed line represents the 
locus of points in the M 2 vs. M x plane, for which the measured lo q and oo are 
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consistent with being due to GW emission only, i.e. if spin-orbit coupling was 
absent (a — 1). This corresponds to a chirp mass M. ~ 0.3 M Q . 
Inserting the measured quantities in eq. (11) and assuming a reference value 
of h = 3 x 10 50 g cm 2 , we obtain: 

1 - aoo " 1.3 x 1034 - 8 X 10 d200 ' (15) 

In principle, the source may be in any regime, but our aim is to check whether 
it can be in steady-state, as to avoid the short timescale problem mentioned in 
2. Indeed, the short orbital period strongly suggests it may have reached the 
asymptotic regime (cfr. [8]). If we assume a = a^, eq. (15) implies (1 — c^oo) — 
8 x 10~ 3 . 

Once UIM and spin-orbit coupling are introduced, the locus of allowed points 
in the M2 vs. Mi plane is somewhat sensitive to the exact value of a: the solid 
curve of Fig. 10 was obtained, from eq. (14), for a — = 0.992. 
From this we conclude that, if RXJ0806. 3+1527 is interpreted as being in the 
UIM steady-state, Mi must be smaller than 1.1 Mq in order for the secondary 
not to fill its Roche lobe, thus avoiding mass transfer. From (1— a^) = 8 x 10~ 3 
and from eq. (4), k ~ 7.7 x 10 45 (c.g.s.): from this, component masses and 
primary magnetic moment can be constrained. Indeed, k = k(ni, Mi, q; a) (eq. 
4) and a further constraint derives from the fact that Mi and q must lie along 
the solid curve of Fig. 10. Given the value of a, hi is obtained for each point 
along the solid curve. We assume an electrical conductivity of ? = 3 x 10 13 
e.s.u. [59, 8]. 

The values of fii obtained in this way, and the corresponding field at the 
primary's surface, are plotted in Fig. 11, from which \xi ~ a few xlO 30 G cm 3 
results, somewhat sensitive to the primary mass. 

We note further that, along the solid curve of Fig. 10, the chirp mass is 
slightly variable, being: X ~ (3.4 4.5) x 10 54 g 5 ^ 3 , which implies M ~ 
(0.26 + 0.31) M . More importantly, E &1 ~ (1.1 1.9) x 10 35 erg s" 1 and, 

since lU/(l-aoo) = E ( ° rb) ~ 1.25xl0 34 erg s^ 1 , we have E gI ~ (9-=-15) E ( ° rb) . 
Orbital spin-up is essentially driven by GW alone; indeed, the dashed and solid 
curves are very close in the M2 vs. Mi plane. 

Summarizing, the observational properties of RXJ0806. 3+1527 can be well 
interpreted in the UIM framework, assuming it is in steady-state. This requires 
the primary to have /ii ~ 10 30 G cm 3 and a spin just slightly slower than the 
orbital motion (the difference being less than ~ 1%). 

The expected value of /ii can in principle be tested by future observations, 
through studies of polarized emission at optical and/or radio wavelenghts [58]. 



5 RX J1914.4+2456 



As for the case of RXJ0806. 3+1527 , we adopt the values discussed in 1 and 
refer to [8] for a discussion of all the uncertainties on these values and their 
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Fig. 10. M2 vs. Mi plot based on the measured timing properties of 
RXJ0806. 3+1527 . The dashed curve is the locus expected if orbital decay is driven 
by GW alone, with no spin-orbit coupling. The solid line describes the locus ex- 
pected if the system is in a steady-state, with (1 — a) = (1 — a ro ) ~ 8 x 1CP 3 . The 
horizontal dot-dashed line represents the minimum mass for a degenerate secondary 
not to fill its Roche-lobe at an orbital period of 321.5 s. Dotted lines are the loci of 
fixed mass ratio. 

implications for the model. 

Application of the scheme used for RX J0806. 3+1527 to this source is not as 
straightforward. The inferred luminosity of this source seems inconsistent with 
steady-state. With the measured values of 1 uj and w D , the system steady-state 
luminosity should be < 2 x 10 32 erg s _1 (eq. 12). This is hardly consistent 
even with the smallest possible luminosity referred to in 1, unless allowing for 
a large value of (1 — a x > 0.15). 

From eq. (11) a relatively high ratio between the actual asynchronism param- 
eter and its steady-state value appears unavoidable: 

|l-a| ~ 2.2(1 -a^) 1 ' 2 (16) 

The case for a. > 1 

The low rate of orbital shrinking measured for this source and its relatively 
high X-ray luminosity put interesting constraints on the primary spin. Indeed, 



1 again assuming I\ = 3 x 10 50 g cm 2 
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Fig. 11. The value of the primary magnetic moment /xi, and the corresponding 
surface B-field, as a function of the primary mass Mi, for (1 — a) — (1 — a™) = 
8 x 1CT 3 . 



a high value of Nl is associated to W ~ 10 33 erg s _1 . 

If a < 1, this torque sums to the GW torque: the resulting orbital evolution 
would thus be faster than if driven by GW alone. In fact, for a < 1, the small- 
est possible value of Nl obtains with a = 0, from which iv[ min ^ = 9 x 10 34 
erg. This implies an absolute minimum to the rate of orbital shrinking (eq. 6), 
3 _/v| ml "- ) /I a , so close to the measured one that unplausibly small component 
masses would be required for E SI to be negligible. We conclude that a < 1 is 
essentially ruled out in the UIM discussed here. 

If a > 1 the primary spin is faster than the orbital motion and the situation 
is different. Spin-orbit coupling has an opposite sign with respect to the GW 
torque. The small torque on the orbit implied by the measured oj could result 
from two larger torques of opposite signs partially cancelling each other. 
This point has been overlooked by [23] who estimated the GW luminosity of 
the source from its measured timing parameters and, based on this estimate, 
claimed the failure of the UIM for RX J1914. 4+2456 . In discussing this and 
other misinterpretations of the UIM in the literature, [8] show that the argu- 
ment by [23] actually leads to our same conclusion: in the UIM framework, 
the orbital evolution of this source must be affected significantly by spin-orbit 
coupling, being slowed down by the transfer of angular momentum and en- 
ergy from the primary spin to the orbit. The source GW luminosity must 
accordingly be larger than indicated by its timing parameters. 
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5.1 Constraining the asynchronous system 

Given that the source is not compatible with steady-state, we constrain system 
parameters in order to match the measured values of W, oo and cj Q and meet 
the requirement that the resulting state has a sufficiently long lifetime. 
Since system parameters cannot all be determined uniquely, we adopt the 
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Fig. 12. M2 vs. Mi plot based on measured timing properties of RX J1914. 4+2456 
. The dot-dashed line corresponds to the minimum mass for a degenerate secondary 
not to fill its Roche-lobe. The dashed curve represents the locus expected if orbital 
decay was driven by GW alone, with no spin-orbit coupling. This curve is consistent 
with a detached system only for extremely low masses. The solid lines describe the 
loci expected if spin-orbit coupling is present (the secondary spin being always tidally 
locked) and gives a negative contribution to cj . The four curves are obtained for 
W — 10 33 erg s _1 and four different values of a — 1.025, 1.05, 1.075, 1.1, respectively, 
from top to bottom, as reported in the plot. 

following scheme: given a value of a cq. (14) allows to determine, for each 
value of Mi , the corresponding value of M2 (or q) that is compatible with the 
measured W, uj and lo . This yields the solid curves of Fig. 12. 
As these curves show, the larger is a and the smaller the upward shift of the 
corresponding locus. This is not surprising, since these curves are obtained at 
fixed luminosity W and uj . Recalling that (1/a) gives the efficiency of energy 
transfer in systems with a > 1 (cfr. 3.1), a higher a at a given luminosity 
implies that less energy is being transferred to the orbit. Accordingly, GWs 
need being emitted at a smaller rate to match the measured iu . 
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The values of a in Fig. 12 were chosen arbitrarily and are just illustrative: note 
that the resulting curves are similar to those obtained for RX J0806. 3+1527 
. Given a, one can also estimate k from the definiton of W (eq. 4). The 
information given by the curves of Fig. 12 determines all quantities contained 
in k, apart from n\. Therefore, proceeding as in the previous section, we can 
determine the value of m along each of the four curves of Fig. 12. As in the 
case of RXJ0806. 3+1527 , derived values are in the ~ 10 30 G cm 3 range. 
Plots and discussion of these results are reported by [8] . 

We finally note that the curves of Fig. 12 define the value of X for each 
(Mi,M 2 ), from which the system GW luminosity i? gr can be calculated and 
its ratio to spin-orbit coupling. According to the above curves, the expected 
GW luminosity of this source is in the range (4.6 -j- 1.4) x 10 34 erg s _1 . The 
corresponding ratios E SI j E < f rb " > are 1.15,1.21,1.29 and 1.4, respectively, for 
a = 1.025,1.05,1.075 and 1.1. 

Since the system cannot be in steady-state a strong question on the duration 
of this transient phase arises. The synchronization timescale r a = a/ a can 
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Fig. 13. The evolution timescale r a as a function of the primary mass for the same 
values of a used previously, reported on the curves. Given the luminosity W ~ 10 33 
erg s _1 and a value of a, r a is calculated as a function of Mi. 

be estimated combining eq. (5) and (6). With the measured values of 14 7 , uj q 
and w , r a can be calculated as a function of I x and, thus, of Mi, given a 
particular value of a. Fig. 13 shows results obtained for the same four values 
of a assumed previously. The resulting timescales range from a few xlO 4 
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yrs to a few x 10 5 yrs, tens to hundreds times longer than previously obtained 
and compatible with constraints derived from the expected population of such 
objects in the Galaxy. Reference [8] discuss this point and its implications in 
more detail. 

6 Conclusions 

The observational properties of the two DDBs with the shortest orbital period 
known to date have been discussed in relation with their physical nature. 
The Unipolar Inductor Model and its coupling to GW emission have been in- 
troduced to explain a number of puzzling features that these two sources have 
in common and that are difficult to reconcile with most, if not all, models of 
mass transfer in such systems. 

Emphasis was put on the relevant new physical features that characterize 
the model. In particular, the role of spin-orbit coupling through the Lorentz 
torque and the role of GW emission in keeping the electric interaction active 
at all times has been thoroughly discussed in all their implications. It has 
been shown that the model does work over arbitrarily long timescales. 
Application of the model to both RX J0806.3+1527 and RX J1914.4+2456 ac- 
counts in a natural way for their main observational properties. Constraints 
on physical parameters are derived in order for the model to work, and can 
be verified by future observations. 

It is concluded that the components in these two binaries may be much more 
similar than it may appear from their timing properties and luminosities. The 
significant observational differences could essentially be due to the two sys- 
tems being caught in different evolutionary stages. RX J1914. 4+2456 would 
be in a luminous, transient phase that preceeds its settling into the dimmer 
steady-state, a regime already reached by the shorter period RX J0806. 3+1527 
. Although the more luminous phase is transient, its lifetime can be as long as 
~ 10 5 yrs, one or two orders of magnitude longer than previously estimated. 
The GW luminosity of RXJ1914. 4+2456 could be much larger than previ- 
ously expected since its orbital evolution could be largely slowed down by an 
additional torque, apart from GW emission. 

Finally, we stress that further developements and refinements of the model are 
required to address more specific observational issues and to assess the con- 
sequences that this new scenario might have on evolutionary scenarios and 
population synthesis models. 
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